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Abstract 



A soft photon approximation is used to calculate the rates of lepton pair produc- 

Q>^ I tion through virtual bremsstrahlung from both pions and quarks. Standard assump- 

■ tions about the evolution of a nuclear system under collision allow pion and quark 

^3 ■ driven total production to be calculated. Comparisons are made with Dalitz decay 

. of light mesons. These mechanisms are expected to be significant contributors to the 

\ soft dilepton mass spectra one might observe in heavy ion collisions at RHIC and 

' LHC energies. 

ix. 

> 

^ \ Thermal photons and dileptons were among the first tools proposed for the observation 

^ \ of thermalized nuclear matter [^. Leptons and photons have small cross sections with the 
nuclear medium, so they will escape without much further interaction from a region of dense 
hadronic matter. This allows them, at least in principle, to carry information about the 
temperature of the primordial state, which may be a thermalized hadronic gas or a quark- 
gluon plasma (QGP). We shall consider here dilepton production from each, concentrating in 
particular on electron-positron pairs. The dilepton invariant mass naturally subdivides into 
four distinct domains. For masses M above the J/ip peak, the spectrum is expected to be 
dominated by Drell-Yan production P . The spectrum around 1 GeV is dominated by p, uj 
and (f) decays. Dileptons with invariant masses in the 1.5-3.0 GeV region seem to be mostly 
thermal [^. Competing sources in this mass region are DD decays @], Drell-Yan processes, 
and pre-equilibrium emission [^Q. The last region of invariant mass (M < nip) is a region 
of soft lepton production [Q. It is expected that below the p mass a non-perturbative QGP- 
induced dilepton radiation may appear. There are many sources of dileptons with invariant 
masses less than rrip. For instance, there are vr" and rj Dalitz decays, tt^tt^ annihilation 
from a thermalized pion gas, virtual bremsstrahlung with subsequent decay from pions and 
quarks, and finally, quark-antiquark annihilation in a QGP. Recent calculations have 
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shown that for masses near 2m^, the contributions from vr*^ and t] Dahtz decays are two 
or three orders of magnitude above the qq annihilation spectrum. However, it has been 
shown for zero pt dileptons having small invariant masses, that the perturbative corrections 
can be several orders of magnitude larger than the Born term [||. So, we consider this qq 
annihilation spectrum ClS db lower limit only. 

The contribution to the invariant mass spectrum from pion scattering with virtual 
bremsstrahlung shown in Ref. [0] is the same order as the 7r° and t] Dalitz decays in the 
region M < 2m^. However, the authors considered only the processes tt^tt" —>■ n^ir^'j*. Our 
purpose here is to supplement their results by generalizing the formulae to include other 
pion-pion reactions, such as tt^tt^ 7r^7r^7*. Contributions from tt^tt^ —>■ 7r^7r^7* are 
ignored because one expects significant cancellation from like-charge interference §]. Note 
the important difference in the ordering of the charges. Furthermore, we use these general- 
izations along with lattice gauge theory data on screening to estimate pair production from 
quark-quark and quark-gluon scattering processes within a soft photon approximation [ITU 



We imagine the following picture for e~^e~ production through virtual bremsstrahlung 
in thermalized hadronic matter. If the initial temperature of the matter is Tj > Tc = 200 
MeV, then virtual bremsstrahlung from quarks occurs. As the system evolves in time it 
cools. If the QGP-hadron phase transition is first order at T = the mixed phase may be 
realized. In this state there will be dilepton emission from both quarks and pions. Then, for 
temperatures lower than Tc, say Tc > T > Tj ~ 140 MeV, one need only consider radiation 
from pions. Let us now estimate the rates of e^e~ production in the quark and pion phases 
in the pursuit of a time-evolution integrated pair production mass spectrum. 

In the reaction tttt — vrvr one expects radiation of real and virtual photons as depicted 
in Fig. with the external lines being pions and the shaded region representing the strong 



interaction. In the soft photon approximation |]TT|, to which we restrict our discussion, only 
the initial and final states radiate virtual quanta. In this approximation, the differential 
cross section for e^e~ pair production is 



dalJ" I f ^ da. 

' \e-J\ 



dM^ Stt^M'^J ' ' dt 

d^pj^d^p 
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5{M^-{p,+p.fy-^:y^dt (0.1) 



where t is the four-momentum transfer in the vrvr collision, M is the invariant mass of the 
e^e~ pair, p+ and p- are the momenta of the positron and electron respectively, and 
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^l' = Ee<5-^e5- J. (0.2) 
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The sum is taken over all photon polarizations, and the current J in the reaction TiaT^b T^c^d 
is 

J'' = -Qa^ - + + Qd^, (0.3) 
Pa-q Pb-q Pc-q Pd-q 

where Qi = Q.,^J\e\ and is the charge of the i th pion. When one of the incoming pions 
is electrically neutral, Eq. ( p.2| ) can be approximated for \t\ < 4m^ by 
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(0.4) 



The square of the photon's energy, ^q, is then replaced by the symmetrized combination 
E{E^ - M^y/^, with E = E+ + = \p+\ + \p^\. Within the same region of vahdity, 
namely, |t| < 4m^, we calculate |e • for all possible charge combinations in iiaTTh T^aT^h 
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^vr" — > 7r^7r° and tt^tt^ tt^tt^) to be 



J\ 
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where 
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3 go V"^^. 
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(0.6) 



The interference function f{s) behaves in the following manner. In the limit ^/s —>■ 2mTr, 
f{s) —2/3 and as ^/s —>■ oo, f{s) —>■ 0. In the small i/i region (^/i — >• 2mT,) the factor 
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|e ■ J|^±^o, but for ^ oo, |e ■ J\l±^^ = 2\e ■ J\l±^o. 



The momentum integration in Eq. 



37r2 M2 



1|) gives 

y/s — 2m^ 



In 



M 



(0.7) 



where (t(s) is the momentum transfer weighted cross section. If da/dt is a symmetric function 
of t and u, then 



cr(s) = 2a^i(s 



Ami 



- 1 



(0.8) 



The elastic pion-pion cross section is parametrized the same as in Ref. 0. For ^/s < 0.6 GeV 
the chiral model expression is adopted, the p-resonance form of aei is used for 0.6 < a/s < 1.5 
GeV and finally, a constant value of 5 mb is taken for ^/s > 1.5 GeV. The e~^e~ production 



rate through virtual bremsstrahlung from pions is |TT 

2^6 



dN 



d^xdM^ (27r)4 
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■^min/'^ 

Ami 



dM^ 



(0.9) 
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where 2;min = '2m.^ + M, d'^x is an element of four-volume of the hadronic matter, and /Ci is 
the modified Bessel function. In Fig. § we show the rate of e^e~ production from the sum 
of 7r^7r° and tt'^tt^ scattering at a temperature T = 200 MeV. Also in Fig. ^ we compare our 
calculations with the rate of production using the chiral model with quartic pion interactions 
Ci = A(7f ■ 7?)^/4 (A = 1.4 from vrvr scattering lengths) ||T^. Our results reproduce the rate 
from the chiral model in the low mass region M < 100 MeV only if ^/s is restricted to the 
non-resonance region, say y/s < 0.48 GeV. Electron-positron pairs with invariant masses 
larger than ~ 100 MeV are attributed to a/s in the resonance region and beyond. 

The soft photon approximation has the advantage that the electromagnetic and strong 
interaction components of the matrix elements disentangle leaving separate multiplicative 
factors (Low's theorem). We utilize this property again to simplify the contributions to 
the e~^e~ spectrum from quark-quark and quark-gluon scattering. The same collection of 
diagrams from Fig. |l] are considered with the charged external lines being u or d quarks (or 
antiquarks) and neutral external lines being gluons. In total we include six quark-quark (or 
antiquark) and four quark (or antiquark)-gluon diagrams. The electromagnetic amplitude 
is obtained from Eq. ( p.5| ) by replacing the pion mass with the quark mass and by setting 
Qa and to the appropriate quark electric charges. The strong interaction cross sections 
dcqg/dt and daqq/dt are well known in the perturbative vacuum at the one-gluon-exchange 



level to be |jT3| 



dt t 
where 




Cab = ; . (0.11) 



For hot hadronic matter this is clearly inadequate. At finite temperatures, static color- 
electric fields are shielded by quarks and gluons, and color-magnetic fields are shielded by 
non-perturbative effects [l^. To include such effects in a finite-temperature evaluation of 



diagrams like the one in Fig. |I|, we use the color- electric mass m\ = 67ras(t)T^, and the 
lattice gauge result for the color-magnetic mass m\j = 25a^(i:)T^ [jl5[ of the exchanged gluon. 
Then for massless quarks whose three-momenta are initially anticoUinear in the frame fixed 
by the finite-temperature medium, we have 



daab _ ^ ag(^)(2t - m| - m|^)^ 

dt " '^"2 it-miy{t-miy ■ ^""^ ^ 

We check the validity of using massless quarks by assuming aE(t) ~ aMit) and introducing 
finite quark masses. The dilepton production rates change very little, so to this level of 
accuracy, Eq. (|0.12|) is acceptable. Note that in the limit mE,mM — >■ (free space) the 
resulting cross section reproduces Eq. ( |0.10| ). Realistically, is also temperature dependent 



owing to many-body effects. Since we make small |t| approximations throughout, we ignore 
the momentum dependence and adopt the lattice gauge theory renormalization group result 
for the temperature dependence |]16| 



"^^^^ = (33-2Ar,)Mar/Tj' ^^'^^^ 
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with a ~ 8. 

Having the strong interaction cross sections, we can use Eq. ( p.9|) to calculate the rate 
of e~^e~ production through virtual bremsstrahlung in qq qq scattering by replacing the 
pion mass with the quark mass and multiplying by the appropriate spin degeneracy factor 
4. For qg qg scattering we must use 

dN 



d^xdM"^ 




(0.14) 

where this time 2;niin = + M and the spin degeneracy to multiply Eq. ( p.l4| ) is 6. Kine- 
matics also changes the argument of the logarithm in Eq. ( p.7| ) to (^i — mq)/M. 

A choice must be made for the quark mass to be used in Eqs. ( |0.5| - |0-9| ) and ( p.l4| ). Chiral 
invariance is restored for temperatures near the phase transition, so constituent quarks 
become current quarks with masses ruq ~ 5 MeV. On the another hand, in a QGP the 
fermion mass induced by many-body effects is rUg = (27ias/<iY^'^T, which is much 
greater than 5 MeV. Since the quark-mass dependence in the rate is roughly l/m^ due to 
Eq. ( p.8| ) and the temperature dependence is roughly as in Eq. ( p.9| ), an upper limit on 
the rate is obtained using nig = 5 MeV at T = 300 MeV. Similarly, a lower limit on the 
rate is obtained by using = 300 MeV at the phase transition temperature. We show 
in Fig. ^ these lower and upper limits on the quark-driven rate, the rate obtained using 
the induced quark mass at two different values of temperature, and finally, the pion-driven 
contribution at Tc. From this one may conclude that the pion gas is slightly more luminous 
than our medium quark gas at Tc, for the range of invariant masses we consider. The rate 
of production from the current quarks is much larger than the pion rate. 

Total e~^e~ production from competing quark and pion processes is obtained only after 
integrating the rate from Eq. (p.9|) over the space-time history of the nuclear system. Using 



the standard 1 + 1 dimensional expansion in a Bjorken picture [T^JT^, we reproduce the 
7r^7r° results of Ref. [0] from mixed plus cooling phases using their values for the relevant 
temperatures. By including the tt^h'^ reactions we find a 20% enhancement over their 
result. The enhancement is smaller than one would expect, but in our case the cross section 
is (Tei ~ a^^'^ -\- a'^ + a'^^'^ with an equal weighting 1/3 of each channel, which is 
different from Ref. ^ where aei ~ a'^^'" +cr'^ and therefore, the weighting is 1/2 for both 
channels. We go on to choose temperatures Tj = 300 MeV, = 200 MeV, and Tf = 140 
MeV in order to investigate other conditions likely to be found in heavy ion collisions at 
RHIC and LHC energies. In Fig. ^ we present the resulting invariant mass spectra through 
virtual bremsstrahlung from quarks and pions separately. The quark result is the sum of 
contributions from the cooling (Tj > T > Tc) plus mixed (T = Tc) phases and the pion result 
is the sum of contributions from the mixed (T = Tc) plus cooling (Tc > T > Tf) phases. The 
pion cooling-phase contribution is the larger contributor to its final spectrum; whereas the 
mixed-phase dominates the resulting quark spectrum. For a more complete comparison with 
other competing processes, we show also the results from 7r° and rj Dalitz decay as well as 
gg — > 7* — > e~^e~ annihilation in the Born approximation. Using our values of Tj, Tc and Tf, 
the medium quark contribution (from cooling plus mixed phases) and the constituent quark 
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contribution (from the mixed phase only) are both lower than the vr" Dalitz, the t] Dalitz, 
and well below the pion contribution. On the other hand, the current quark contribution is 
larger than 7r° or t] Dalitz; it is even larger than the resulting pion contribution by a factor 
of 10 or more. However, we believe the true quark contribution will be nearer to the medium 
quark result. This being true, virtual bremsstrahlung from pions will be the largest source 
of low-mass e~^e~ pairs. This points out the need for more elaborate calculations, perhaps 
based on nonequilibrium approaches. 

Assuming only soft virtual photons greatly simplified our task by eliminating the need 
to exactly evaluate the many Feynman diagrams that contribute to e~^e~ production. But 
the validity of this popular approach may be questioned. For instance, the restriction that 
the momentum transfer \t\ be less than 4m^ for pions and less than Am'^ for the quarks 
is arguably too restrictive. Also, many-body effects introduce temperature dependences 
into the meson masses, widths and coupling constants and therefore, into the tttt elastic 
cross section likely to be quite strong near T^. Detailed discussions of these effects will be 
presented in a future paper [|19| . 
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FIGURES 



FIG. 1. Lepton pair production through virtual bremsstrahlung in scattering of particles 
a b ^ c d. The particles might be charged or uncharged pions, quarks or antiquarks or gluons. 
The shaded region indicates a strong interaction. 

FIG. 2. Chiral model predictions for the e~^e~ production rate as compared with our calcula- 
tion. The dashed lines are our estimations with varying values for the upper limit on the energy 
and the solid line is from the chiral model. 

FIG. 3. Mass dependence of quark-driven and pion-driven production rates at temperatures 
T = 200 MeV and T = 300 MeV. 

FIG. 4. Dilepton mass spectra for LHC energies: The solid curve is the total contribution 
from pion scattering processes, the dot-dashed curves result from quark-quark and quark-gluon 
scattering with different quark masses, the long-dashed curve is the result of qq annihilation in a 
Born approximation, and finally, the short-dashed curves are Dalitz decay spectra. 
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